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Abstract We analyze constrained quantum systems where the dynamics do not 

Q , preserve the constraints. This is done in particular for the restriction of 

a quantum particle in R n to a curved submanifold, and we propose 

a method of constraining and dynamics adjustment which produces the 

right Hamiltonian on the submanifold when tested on known examples. 

j^j This method will be the germ of a "Dirac algorithm for quantum con- 

r^ ■ straints." We generalise it to the situation where the constraint is a 

general selfadjoint operator with some additional structures. 
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1. Introduction. 

The problem of enforcing constraints in classical mechanics has had satisfactory 
solutions for some time now [Di, GNH, MW]. The same is not true in quantum 
mechanics, contrary to claims of "quantizations" of classical constrained systems 
after constraining. Quantization maps are rarely unique nor well-defined [VH, Wo, 
Ri, GGT]. Moreover, using noncommutativity, some constraints can be defined in 
quantum mechanics with zero classical limit. This means that the general problem 
of quantum constraints needs to be solved in the quantum arena without appealing 
to classical methods. 

There is presently a wide variety of methods for imposing quantum constraints 
[Di, GH, HT, La, SW] but these methods deal with the kinematics only at the 
quantum level, and their interrelations are unclear. When the given quantum dy- 
namics preserves the constrained subsystem that is fine, one only needs to restrict 
it to the constrained subsystem to obtain the constrained dynamics. However, 
when this is not the case the problem arises of how to appropriately adjust the 
constraint set and the dynamics into a stable system. (New constraints produced 
by such a method are generally called secondary constraints). The obvious strat- 
egy of fixing the dynamics and extending the constraint set to its orbit under the 
dynamics, produces the wrong physics in examples. It is our opinion that the 
question of secondary quantum constraints and dynamics adjustment is a physical 
problem, and cannot be decided by mathematical arguments alone. In this paper 
we set out to solve this problem of secondary quantum constraints in the limited 
class of systems consisting of the restriction of a quantum particle on R n to a 
smooth (possibly curved) submanifold i.e. the case of holonomic constraints (cf. 
p75 [Ar]). The motivation for this choice is as follows: 

(1) there are several examples of quantum mechanics on surfaces available (e.g. 
S 2 ) against which we can test the results of our analysis. 

(2) in R n there is ample geometry available to guide the intuition, e.g. a metric, 
hence normal vectors and Lebesgue measure, 

(3) generalisations are easy from this class, in fact in Sect. 6 we propose a gener- 
alisation of it to the case of a general selfadjoint operator on a Hilbert space 
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with additional structures, 

(4) a constraint in this situation is of the form (p(q) = which has an immediate 
obvious quantization in the Schrodinger representation, whereas a constraint 
of the form cp(q, p) = need not have a unique quantization, 

(5) physically these systems in IR 3 can be approximated by very thin films or 
wires and so are close to experimental verification. 

(6) the exposition of constraint methods for these systems are more transparent 
than for more general ones, 

(7) holonomic quantum constraints have been considered in the physics literature 
[Ma, Fu]. 

The kinematics of quantum systems of this type has already been done concretely 
by Landsman [La] and abstractly at the C*-level by Grundling and Hurst [GH], so 
any new proposed method should be compatible with these. We will not consider 
BRST-methods [HT], as we believe these not to be equivalent at the quantum 
level to the other methods [Gr, LL, Mc]. So to summarize; we aim to solve for 
this class of systems the problem of secondary quantum constraints, i.e. given a 
constrained system where the dynamics do not preserve the constraints, to find a 
method of adjustment which agrees on known constrained systems with the right 
physics. 

The architecture of the paper is as follows; in Sect. 2 we summarize the 
Dirac procedure for constraining a classical particle to a submanifold, both local 
and global, and we also do two examples. In Sect. 3 we consider the question 
of quantum constraints, we first summarize the usual Dirac procedure, then start 
to analyze the situation of restricting a quantum particle in M. n to a lower 
dimensional submanifold F . We obtain the constraining map k from the 
Hilbert space L 2 (M n ) of the unconstrained particle to that of the constrained 
particle L 2 (F) , and we discuss generalisations to other types of constrained 
systems. In Sect. 4 we solve the problem of how to obtain from the Hamiltonian 
of the original particle a constrained Hamiltonian on L 2 (F) , and we apply it to 
two examples;- restricting a free particle in M 3 to a sphere and a cylinder, and we 
obtain in each case the correct Hamiltonian. In Sect. 5, we consider how to obtain 
observables on L 2 (F) from the original observables on L 2 (M. n ) , i.e. we analyse 
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when we obtain sensible liftings of operators through the map k , and how to 
obtain a suitable field algebra for the constrained particle. In Sect. 6 we suggest a 
generalisation of the method of the preceding sections to the constraining situation 
where the constraint is a selfadjoint operator on an abstract Hilbert space, with 
zero in the continuous part of the spectrum. We find that we need to assume 
some additional structure to do the work which the geometry of M. n does in the 
method of the preceding sections. The reader in a hurry can start with Sect. 3. 

2. Classical Secondary Constraints. 

We start by recalling the basic Dirac-Bergman method [Di, Su, SM] which is of 
course a local procedure. However to keep this brief, we only present it in the 
context of a particle in M. n constrained to an (n— 1) -dimensional submanifold 
r , which is the system at the focus of this paper. 

The full phase space is M. 2n with generic point (q, p) and the usual Pois- 
son algebra V = (C°°(IR 2n ), {•,•}) is the setting for the problem. We specify 
a smooth (n — 1) -dimensional submanifold F in M. n as the zero set of a 
smooth function ip G C°°(IR n ) , i.e. F = </? -1 (0) . We also assume Wip ^ on 
r so we can use the gradient to define a normal on F . The Hamiltonian for a 
particle in a potential V is: 

2m 
So we have a constrained system with a single primary constraint y(q, p) := <£>(q) 
and primary constraint manifold in phase space </? -1 (0) =Tx R n =: F p . The 
Dirac procedure starts with the new Hamiltonian: 

h p '■= w~ p • p + ^(q) + W> > 

where \x G V is a free multiplier, to be determined below. So the new time 
evolution for A G V on F p is 

A \ F p = {A, H p } r T p = ({A, H c } + n ■ {A, ^ \F P . 

Since time evolution must preserve F p we require that 

= ip \F p = {p, H p } \F p = ^p • ( Vv?) \F p . 
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Thus x := ~P ' (Vy) is a secondary constraint determining a smaller submani- 
fold Tf := </? _1 (0) fl x _1 (0) , which again must be preserved by time evolution. 
In the examples in the rest of this paper we will constrain free motion to a surface, 
so for this case (V = 0) : 

i\T f = = ({£, H c } + [i ■ {£, y?}J [T/ , where £ denotes x, 99 

so = ({i P • (v v ), ^p • p} + p{£p • (Vv?), *>}) rr y 

= (^p • (Vv, p • p} + ^(Vv?) • {p, ¥>}) rr y 
= ( Atp • (v( P • VyO) - £/*(Vv>) • (v*>)) rr, . 

Since V</? 7^ on T , we solve on Tf : 



fi 



p-(v(p-vy)) _ (p-v)V 

m\\V(p\\ 2 m||Vv?|| 2 

(P-V)V 



thUS ^=2^P-P + ^ m || Vy ,|| 2 • 

Notice that the second term in ii/p when extended near T/ acts like a "poten- 
tial" to keep the motion on F , and that the secondary constraint 

x=(i P -v^r = o (2.2) 

guarantees that the motion is tangent to the surface and there are no more con- 
straints. Moreover, since 

{ip, x}rr / = ^||v v || 2 rr / #o, 

the constraints are second class. The Dirac algorithm continues to construct a 
Dirac bracket, but we will not need this. For later use we specialize the above to 
a sphere and a cylinder. 

Let r be a sphere of radius a in IR 3 . The phase space is R 6 , T p = 
S 2 x M 3 and we have 

v?(q) = q • q - a 2 , so x = ^P • Vy? = ^p • q 

and since V(p • Vy?) = 2p here, we get 

1 (, . q-q-a 2 

— p-p(l + — jj— 
2?n V ||q 



H p = —P -P(l+ ,^.,,2 ( 2 - 3 ) 
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An explicit calculation produces 



where L = q x p 



H\T f = ^-p-p i = -^—^ (2.4) 

pi / 2mf f |„|=a 2ma 2 v ; 

qp = 



is angular momentum on the sphere. 



If we take F to be a cylinder of radius a around the z-axis in R 3 , we 
have 

^(q) = Qi + Q2 ~ o 2 , so X = ^ (PiQi + P2q2) 

ffp = ^ p . p + (gf + gi _ o2) _|+^_. (2 . 5) 



And now: 



1 /£ 



2 



^Ir/ = AP-P|^_. =2^(i+P3)tr/ (2.6) 

P191=-P2<32 

where L 3 := Qip 2 — ?2Pi is the usual angular momentum around the (/3-axis . 



Next we recall the global version of the Dirac-Bergman method of constraints 
as worked out by Gotay, Nester and Hinds [GNH], which is a considerable sim- 
plification. Start with a symplectic manifold (M, a;) (finite dimensional here) 
and a Hamiltonian H G C°°(M) . Let M\ C M be a given primary constraint 
submanifold. Now the evolution equation i x uj = dH may not have solutions for 
the vector field X tangential to Mi , so let 

M 2 := { m e Mi | dH{m) = i x u{m) for some X e T m M 1 } . 

However the time evolution must now preserve M 2 , so iterate: 

Mfc_|_i := {m E Mk | dH(m) = i x u(m) for some X E T m M^ } . 

When this iteration converges sensibly, we obtain a final constraint manifold M c 
on which dH = i x uj has solutions X G r°°(M c ) = smooth vector fields on 
M c , which provide the desired time evolutions. Geometrically we can think of 
M c as the largest submanifold which has a vector field X which is tangential at 
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each m G M c to a trajectory of the time evolution of the original unconstrained 
manifold M . The completeness of X is still an open question for the general 
analysis. (Below we will refer to the above as the GNH-algorithm). 

We will take this geometric picture as the guiding principle for constraining 
the dynamics of a quantum system below in Sect. 4. 

3. Constraining the representation space. 

For reference we recall the usual Dirac prescription for first class quantum con- 
straints. Given a field algebra A preserving a dense domain A in a Hilbert 
space 7i , as well as a Hamiltonian H G A sa and a set of constraints 
{^e A sa I i G I } (where I is an abstract index set) , one selects the physical 
subspace by 



H phys :={ije A | &V = Vi} 

and this is assumed to be nonzero (hence zero must be in the discrete spectrum 
of each (pi ). Decompose Ti = 7Y p h ys © ^ph ys > so ^ or ea ch A G A we have 
a decomposition A = I , J . Clearly the only A G A which, together 

with A* can be restricted to 7i p h ys are those of the form A = I _, 
Denote the *-algebra of these by O . Those which restrict trivially are of the 
form A = I n , J and comprise an ideal T> of O . Then the *-algebra of 
physical observables are 

o\n phys ^o/v=-.n. 

If exp itH preserves H p h ys then Ad exp(itH) lifts to automorphisms on 
1Z . Otherwise we need to alter it (cf. Sect. 4). Let -Pphys be the projection 
on 7iph ys and define the sesquilinear form 

C01, tp2)D ■= (PphysV'l, ^PphysV^) = (lj>l, Pphys V^) (3-1) 

then Ker (•, -)d = ^ph y s an< ^ we wr he -P p h ys as a factorization -P p h ys : 7~L — > 
Tt/Ker (■, -)d — ^ P hys (obvious identification), and we can think of O as all 
A £ A which together with A* can lift through the factorisation. 

Next we introduce our main object of study. Consider a quantum particle in 
R n in the Schrodinger representation. Its basic data consists of the Hilbert space 



H. = L 2 (IR n , fj) (with (i the Lebesgue measure), and on the dense subspace 
C^°(W n ) of smooth functions of compact support, we have the position operator 
(qi/>)(x) := xt/?(x) , the momentum operator (p)(x) := (zVV0( x ) an d a Hamil- 
tonian H (which is ^p ■ p when the particle is free). These operators are all 
essentially selfadjoint and preserve C^°(M. n ) . We wish to constrain this particle 
to a given smooth (n — 1) -dimensional submanifold rcR" . Assume we have 
a bounded real-valued constraint function ip G C°°(IR n ) with Wip nonzero and 
bounded on a neighbourhood of F , and 

r={xGl n | (p(x) = 0} =^ _1 (0) . 

(Notice that the level hypersurfaces of <p form a foliation of this neighbourhood 
[To].) Quantize <p by the multiplication operator y?(q) , i.e. 

(^V)(x) = v?(x)-V(x) V^eH 

which is selfadjoint since <p is real-valued, and bounded because <p is. 
Assume we are given a concrete field algebra T C B(Ti) , which should be a C*- 
algebra containing all relevant operators in bounded form, e.g. ip , exp za • q , 
exp ia ■ p ( a G M. n ) . The choice of field algebra will turn out to be important 
for producing a meaningful algebra of observables in the constrained system, but 
we will return to this matter in a later section. In summary, we are assuming the 
following data: 

• the operators p , q and H on C^°(K n ) C L 2 (R n ) , 

• Abounded ip G C°°(IR n ) with r = </? -1 (0) a smooth (n — 1) -dimensional 
submanifold, and with Wip nonzero and bounded in a neighbourhood of F , 

• A unital field algebra T on L 2 (IR n ) . 

We would like to find a way of imposing the constraint " ip = "on this system. 
However since F is a null-set of // , we have Ker ip = {0} , and so Dirac's 
method of restriction to the kernel of the constraint fails. The fact that zero is in 
the continuous spectrum of <p is physically due to the uncertainty principle;- a 
particle confined to F will violate Heisenberg's uncertainty principle in a direc- 
tion normal to F . 
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At present there are two methods which can handle the above situation;- 
the T-procedure of [GH], as well as Landsman's application of Rieffel induction 
[La] . The T-procedure "ignores" the original representation of JF , considers T 
as an abstract C*-algebra, for which it then finds those representations tt for 
which zero is in the discrete part of the spectrum of Tc(ip) and works out the 
algebraic structures associated with factoring Ker n((p) out of each of those rep- 
resentations. On the other hand, Landsman's method assumes a locally compact 
group H with a continuous proper action on M. n such that T is precisely 
the set of points left invariant. Assuming some initial representation n of H 
on L 2 (M. n ) , he makes C c (lR n ) into an (A— £>)-bimodule where B is the 
algebra C C (H) (with convolution for multiplication) and A is the part of the 
commutant of tt(H) preserving C c (M. n ) . There is then a natural rigging map, 
and this allows a Rieffel induction from the trivial representation of B to A . 

The T-procedure selects the set of all representations in which the Dirac al- 
gorithm makes sense, which necessarily excludes the original representation of the 
current system. This means that one loses the original physical interpretation of 
the state vectors, and one has the problem of which representation to choose for 
the constrained system. From the point of view of the physics, one may some- 
times be more interested in constructing the constrained operators directly out 
of the original ones. Landsman's method does this, but requires some additional 
group structure, and quantizes a fairly small algebra. Both methods do only the 
kinematics, and if the dynamics do not preserve T it is hard to see from these 
structures what should be done. 

We will propose a method for constraining the kinematics of a quantum par- 
ticle to r which generalises the usual Dirac prescription and in which the con- 
strained operators are explicitly constructed from the original ones. It will then 
be fairly easy to see what to do with the dynamics. 

In a sense the solution of the kinematics problem is already known. One takes 
the subspace C c (M. n ) C L 2 (M, n ) and restricts it to T , thus obtaining a dense 
subset of the physical Hilbert space L 2 (T, 7) ( 7 denotes the measure induced 
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on r by p ). The restriction is the same as factoring out the subspace 

N := { f G C c (R n ) I f\T = 0} 

from C c (W l ) . The observables are those operators which "restrict to F ," i.e.. 
operators A which preserve C c (IR n ) and N , in which case A lifts through 
the restriction to define an operator on C C (T) C L 2 (T) . This is the procedure 
which we intend to refine. Note however that the restriction map p : C c (IR n ) — > 
C c (r) is unbounded with respect to the L 2 -norms and it is not closable as an 
operator (let ipk £ L 2 (R n ) be a sequence of continuous functions with compact 
support, all of which have the same restriction to a nonzero 9 G L 2 (T) , and 
which converge to zero w.r.t. the L 2 -norm . Then ip^ — »■ , p(iftk) ~~ *■ $ 7^ , 
so \im p(ifjk) t^ 1 /o(lim^fc) ). So we expect some pathologies to arise in operator 
questions. We would also like to build the inner product of L 2 (T) explicitly out 
of the inner product of L 2 (R n ) . To do this, we adapt a well-known argument 
from statistical mechanics [Kh], and p80 of [Fa]. 

Since V</?(x) ^ on a neighbourhood of F , the critical points of ip are 
a finite distance away from T . So we have the geometrically obvious: 

Lemma 3.2. There is a t > such that in </? _1 [— t, t] we have locally a 
smooth curvilinear orthogonal coordinate system (y, ip) where 
y = constant are curves normal to the level hypersurfaces of 

V? . 

Proof: Consider the vector field X(x) = n(x)/|V</?(x)| on tp~ x (— t, t) , where 
n = V</?/|Vy?| . Then its integral curves are normal to all level surfaces 
of ip in f~ l (— t, t) , and moreover for each s G (— t, t) we have 
a diffeomorphism a s : </? _1 (0) — > </? _1 (— t, t) such that ao = t and 
^a s (x) = X(a s (x)) [AM]. In fact p(a a (x)) = s for all x G </? _1 (0) = 
r , which we see as follows. For small e : 

/ x n ( X ) ^/ \ 

so V?(a £ (x)) = (p(pc) + (e/|Vp(x)|) V<^(x) • n(x) + 0(e) 

= V?(x)+e + 0(e) VxG <^ _1 (— t, t) 
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using a Taylor expansion and n = Vy/|V<^| . So 4-(p(a s (x.)) = 1 
with solution </?(a s (x)) = s + </?(x) and if x G V , the last term is 
zero. Thus we have diffeomorphisms a s : </? _1 (0) — > </? -1 (s) . Now equip 
an open neighbourhood [7 C (/? _1 (0) with local coordinates y by a 
chart, then a s will equip a 8 (U) C </? _1 (s) with the same coordinates. 
Thus we have the desired local coordinates (y, <p) on the set V : = 
<f~ 1 (—t, t)C\a ( _ t t JU) . (The incompleteness of X on ip~ 1 [—t, t] is 
not an obstacle to the argument;- just attenuate X with appropriate 
smooth bump functions to make it integrable with speed 1 on ip where 
we want it). | 

For a function / G L 1 (M n ) we have 

/ /(x)rf^(x)= / f(y,(p)J(y,(p)dydip 
Jv Jv 

where J is the Jacobian and V is as in the preceding proof. Observe that for a 
fixed if = s that J(y, s) dy is not yet the surface measure d"f s on </? _1 (s) . 
For that the orthogonal coordinate ip needs to be expressed in terms of the length 
of the curves y = constant . Since d(p — \W(p\ d£ with d£ the length measure 
on a curve y = constant , we conclude that d^ s (y) = [J ■ |V<^|](y, s) dy . Now 

/ |V^(x)|/(x)^(x)= / /(y,^)[J-|V^|]dyd^ 
Jv Jv 

= J (J f-J-\V<p\dy)d*p (3.3) 

and this is the expression we wish to exploit. Let / G C c (lR n ) , and let the 
thickness of the shell <p~ 1 [—t, t] around F approach zero, then 



lim — 



vrup- 



\Vip\-fdn=hmJ- f ( [ f-J-\Vcp\dy)d<p 
s^o 2s J_ s V Ju J 

= [ {f-J-\V<p\)(y,0)dy 

Ju 



where the use of the fundamental theorem of calculus is justified because the 
function 



'u 



[ (f-J-\V<p\)(y,<p)dy 
Ju 
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is continous due to the uniform continuity of the integrand;- a consequence of 

/ G C c (R n ) and J • \V<p\ G C°°(V) . Thus 

Ihn-L f \V<p\.fd»= f(f\T)dry. 

By doing this for all open patches U C Y (equipped with charts) , we conclude 
liml/ \V<p\-fdfi = f(f\T)d 1 (3.4) 

for all f eC c (R n ) . 

In particular, let z/>i, ^ G C c (M n ) , then we have 

l (r) t— o 2i y^-if-t,*] 

1 

= limi(^,^ 2 ) (3.5) 

where we used the notation (/m/>)(x) :— |V<p(x)|^(x) and P t := x - 1 \-t t] an< ^ 
this is the desired relation between the inner products of the initial space L 2 (R n ) 
and of the physical space L 2 (T) . (Note that P t is just the spectral projection 
of ip on the interval [—t, t] , and that h is a positive bounded operator on 
PtH ). Now the right hand side of this equation will exist for a much larger class 
of functions in L 2 (R n ) than C c (lR n ) , though for some of these the restrictions 
to r on the left hand side may not be defined. However the limit in (3.5) will 
definitely fail to exist for some elements of L 2 (R n ) . Define 

(V>1, V> 2 )r := gm ^ (Vi, hP t fo) (3-6) 

for all pairs ipi, ip2 for which the rhs is defined and finite. To obtain a positive 
sesquilinear form from (3.6), we need to decide on a dense domain on which (•, -)r 
is defined and finite. Clearly C c (M. n ) is one such domain, but it is not maximal. 
In fact, (•, -)rfC c (lR n ) has no closed extensions as a sesquilinear form (cf. [RS] 
p278), so there is no canonical way of getting a maximal domain. The set 

S := { V> G L 2 (R n ) j (V>, V)r exists and is finite } 
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is not a linear space since (ipi, V^r need not exist for ipi, ip2 G S as one can 
verify with easy examples. Nevertheless, some domains containing C c (lR n ) seem 
quite natural, e.g. the domain C^ consisting of those V £ £ 2 (K n ) for which 
there is some t > (depending on ip ) such that Ptip G C c (</? _1 [— t, t]) , and 
clearly the argument which led to eq. (3.5) is still valid for ipi G C v . Now C v 
is still not maximal and we can of course choose a maximal defining subset for 
(•, -)r in S containing C v , but this choice is not unique. Equip the factor 
space C^/Ker (•, -)r with the inner product (•, -)r , making it into a pre-Hilbert 
space which we complete to obtain the Hilbert space Hr • We identify Hr with 
L 2 (r) in a natural way: 

Theorem 3.7. With the structures above, denote by C : C^ — ► C^/Ker (•, -)r 
the factor map, then 

(1) (Ci; 1 ,Cij 2 )r={^i\T,i; 2 \T) for all fa G C c (R n ) , 

L (T) 

(2) there is a unitary U : Hr — ► £ 2 (r, 7) such that C/C'0 = ip |T 
for ah V e C c (M n ) , and (t/C)(C c (M n )) = C c (r) . 

Proof: (1) Since (V>i, "02 )r = (CV'ij C?/>2)r , this is already proven above. 

(2) Let V e C c (M n ) and define t/C^ := if>\T which is in C C (T) 
since T is a closed subset of M. n . To see that U is well defined 
on C c (R n )/Ker (•, -)r , let Vi, V>2 G C c (M n ) with (V>i - ^2) fr = . 
Then by (3.4) we find ^-^G Ker (•, -)r , i.e. Cipi = Cip so U is 
well-defined. Moreover from (1) we see that U is unitary on the dense 
subspace C[C c (M n )] hence it extends to a unitary on Hr • Clearly 
C c (lR n ) |T = C c (r) , so since the image of U contains a dense subspace 
and U is unitary, we find U : Hr — ► £ 2 (r) onto. | 

This means the restriction map can be written as k := U o C : C v — ► £ 2 (r) , 

but now it implicitly involves a limiting process, which we will exploit below. 

Note that C c (R n ) n Ker (•, -)r = { $ G C c (M n ) I Vfr = 0} and that « is 

unbounded (w.r.t. the Hilbert space norms). 
Example. We want to see what the above structure looks like in the restriction 
of a quantum particle in R 3 to a sphere. So take Ti. = L 2 (M. 3 ) , 
pip = iVip , qV'W — x , 0(x) as usual and for the constraint <p(x.) = 
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|x| 2 - a 2 , a > . Now (V<p)(x) = 2x 7^ except if x = , and 
so the critical point of ip is distance a away from F = </? _1 (0) . We 



have if 1 [— s, t] = { x G 



xr G a' 



s, a 



+ *]} 



, s, £ > . So 



(V>i, V>2)r = Jim glCV'i. ^^2) 



d_ 
dt 



v~ 



-t,t] 



2|x|(^ 1 ^ 2 )(x)(i 3 x 



V^+t 



d 



r 2 dr I d£lr (V'iV'2) (3 



t=o 



t=o 



where we used polar coordinates and O denotes the measure on the 
unit sphere S 2 . For functions of the form ipi{x) = fi(r)£i(0, 0) , 
£i G L 2 (S 2 ) , and fi continuous on (a — s, a + e) we have by the 
fundamental theorem of calculus 



(V'l. i>: 



2)T 



d_ 

~dt 



v / a 5 +t 



y/aP^t 



r 3 fif2{r)dr 



t=o 



66(0, 0)^ 



a 2 fi(a)f 2 (a)J^ 2 dn. 



as expected from (3.5). Moreover we have that ifti G Ker (•, -)r when- 
ever fi(a) = . 

One may ask how the method above should be generalised to bigger classes of 

constraints, so towards that a few remarks. 
Remarks. (1) If F has corners (so ip is not smooth) we can develop an "approx- 
imation" to it by a sequence of smooth submanifolds {I\} , but need 
to choose a notion of convergence for Fk — *• T . 
(2) If r has edges we can combine the procedure above with a Dirac 



constraining. For example if we want to constrain a particle in I 
the upper hemisphere r + of the sphere F of the last example, 



to 



+ - = 



{ xe 



|x| = a, x\ > } 



xl 2 — a 2 to 



then first apply the method above to the constraint y(x) 

obtain a dense subspace of L 2 (F) , and follow this by applying a Dirac 

procedure via (3.1) to the constraint £(6>, (/)) = Xi 7r /2i(^) m L 2 (F) • 
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(3) To constrain a particle in MJ 1 to a smooth (n — k) -dimensional 
submanifold V , we need k independent constraints cpi,...,(pk such 
that r = (fi (0) fl . . • Hy?^ (0) with critical points well away from T . 
Now it is clear how to adapt the method based on (3.2);- we choose a 
curvilinear coordinate system (y, <pi, . . . , </?&) where <pi are the coor- 
dinates measured along the trajectories of the gradient fields V</?i , and 
y are the coordinates along the manifolds <p^ (ci) D • • • fl <p^ (ck) , c 
in some e-neighbourhood of zero in M. h . Then with J the Jaco- 
bian, the measure on F is d'j — (J • |V</?i| • • • |V<^fc|)(y, 0, . . . , 0) dy 
so (3.3) becomes for / e C c (R n ) : 

\im[2 k s 1 ---s k }- 1 [ \V<pi\---\V<p k \-fdn= !{f\V)d 1 

s ^° 7^5-1 [-s,s] JT 

where </? _1 [— s, s] := ip^ 1 [—si, si] fl ■ • • fl (p^ 1 [—Sk, Sfc] • Further adap- 
tations are straightforward. 

Alternatively, we can impose the constraints ipi one-by-one, but for 
this we need to generalise the procedure above for M. n to general Rie- 
mannian manifolds. 

(4) For ip which are positive, T = </? -1 (0) consists only of critical 
points, though if other critical points are well away from F , there is 
still a neighbourhood of F foliated by the level hypersurfaces of ip , 
and it is possible to adapt the method above to this situation. 

(5) In the case where <p -1 (0) = F U A with ^(r) =0^ ,u(A) , 
mA = , neither the Dirac method (which will only produce L 2 (A) ) 
nor the method above seem appropriate. So what we want is a method 
which will give the method above on F , and the Dirac procedure 
on A . (We choose not to allow surface terms on A ). Assume 
the critical points of ip are well away from F , and that A is 
the closure of an open set. That is, we want a densely defined map 
k : L 2 (R n ) -► L 2 (F) ©L 2 (A) . Recalling that L 2 (A) = P {0} L 2 (R n ) , 
this can be done by defining: 

(V>i, ^2) R := (Vi, P {0} V 2 ) + Jiml«, hP t il>?) 
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for all ipi £ T)r := { ip G H | ip N £ C_p} where ip N := x^ • V' an d 
N is an open neighbourhood of F such that N fl A = . Then 
proceed as before, taking D^/Ker (•, •)# with k the composition of 
the factor map and the unitary from D^/Ker (•, -)r to L 2 (r)©L 2 (A) . 
(6) At the cost of complicating the current exposition, we can enlarge 
the domain C v of k considerably. For instance, by splitting the limit 
(3.6) into the average of the limits from above and below in t , and 
letting the domain consist of those V £ L 2 (M. n ) for which ip — > ip(y, (p) 
is continuous from above and below at zero (where we used the curvilinear 
coordinates of Lemma 3.2), we obtain another useful domain for k .We 
chose not to do this, because the additional analytic details would have 
obscured the simplicity of our central idea. 

(5) The methods proposed here for dynamics reduction can be extended 
to other holonomic constraints, i.e. if we are given a constraint <£>(q, p) 
which involves only one member of each canonical pair (^, pi) , then 
a partial Fourier transform can convert it to the type y?(q) considered 
here. So there is a unitary transformation of such a system to one of the 
present type. 

4. Constraining the Dynamics. 

In this section we continue the analysis of the problem of constraining a quantum 

particle in M. n to a subset and address the problem of how to constrain the 

dynamics, i.e. how to construct out of the given time evolution on L 2 (M. n ) an 

acceptable time evolution on the physical Hilbert space. We consider four cases: 

(1) For an ordinary Dirac constraining we restrict for example a quantum particle 

in R n with Hamiltonian H to live on a set T which is the closure of 

an open set S (hence /u(T) ^ ). Kinematically, one constrains via the 

projection P p h ys : L 2 (M. n ) — > L 2 (T) , as described at the start of Sect. 3. If 

the dynamics preserve T , i.e. each Ut '— expitH is of the form Ut = 

( jj ° h J with respect to the decomposition L 2 (R n ) = L 2 (T) © L 2 ^) 1 - , 

then there is no problem;- one just restricts Ut to L 2 (T) to obtain the 
constrained dynamics. In practice this is a rare occurrence. 
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(2) For the case of a Dirac constraining as in (1) where the dynamics does not 
preserve T , e.g. if the particle is free, H = p 2 /2m and T is compact, 
then H\L 2 {T) is not selfadjoint. To equip the constrained particle with a 
time evolution, we need to choose some selfadjoint extension of the symmetric 
operator H\C^°(S) (for this to make sense we need to assume that H 
preserves C^°(S) , which will be true if it is a differential operator). This 
amounts to the choice of boundary conditions, i.e. deciding how the particle 
should behave at the walls (e.g. reflection), and this is a physical choice which 
cannot be determined from mathematical considerations alone. 

(3) For the problem of the last section where we constrain a particle in M. n with 
Hamiltonian If to a lower dimensional submanifold F , we will need to 
assume that H is "smooth enough" near T , i.e. H (DomH fl dp) C dp 
and ^(Dom H T)C '<£,) is dense in L 2 (T) . In the case where H "preserves" 
r , i.e. H preserves DomH PlKerK , H will clearly lift through k 
and we only need to define the constrained Hamiltonian Hy on L 2 (T) by 
HrK(ip) := K,(Hip) for all ip G DomHnd,p , and require that it is essentially 
selfadjoint. 

(4) The case of constraining a particle as in (3) but where H does not preserve 
r , will be our object of study for the rest of this section. By analogy with 
case (2), we expect that some physical choices will need to be made. When 
H is a Schrodinger operator with smooth potential, physicists already know 
what to do;- one lets the constrained Hamiltonian be 75— Ar + V IT where 
Ar is the Laplacian on F . 

Henceforth we assume that H on L 2 (M n ) does not preserve V in the sense of 
(3) above. Recalling that k : d v — ► L 2 (T) involves a factorisation, it is natural 
to look for linear sections of k which one can use to construct liftings of operators 
on d v to £ 2 (T) . The image of such a section deserves a name of its own: 
Def. A transverse space is a linear space T C d^ such that Ker k fl T = {0} 
and k(T) is dense in £ 2 (r) . 
Transverse spaces are abundant, and below we will use physical arguments to 
make a choice of a transverse space but first we show how one can use a transverse 
space to constrain H :- given a transverse space T C Dom H fl d^ such that 
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HT C Cp , we have the options: 

(i) Define the constrained Hamiltonian H^ with dense domain k(T) by 

This is well-defined since by definition of transverse spaces, k is injective 
on T . 
(ii) Given that the choice of T has some physical content (see below), one may 
object to (i) in that H is not forced to preserve T in any sense. If 
we want to include such a restriction, since it is only the behaviour near F 
which should be important, we can restrict H to T "in the limit," i.e. we 
consider the projection Plj-H as t goes to zero in the definition of k , 
where Pj- denotes the projection onto PfT . So we propose an alternative 
constrained Hamiltonian H^ on k(T) by 

(H (ii) K(i>), «(0)) L2(r) = Jim i(i4#V>, hP t (j>) 

for all ij), (j) eT . 
Below we will test both methods for a suitable choice of transverse space. This 
choice is the issue we now want to address. Recall that in the classical global 
method of Gotay, Nester and Hinds [GNH], the idea was to adjust the dynamics 
so that the motion is always tangential to F , thus forcing the particle to remain 
on r . We look for some quantum mechanical version of this. 

First we want to give meaning to the concept of motion "tangential to F 
for a quantum particle. Whilst the classical state of a particle is a point in phase 
space, quantum mechanically a state is here a vector ip in L 2 (IR n ) . Classically, 
a particle with constant mass at position q moves tangentially to F if its 
momentum p = mq is tangential to F , i.e. p G T q r . Inspired by this, in 
a quantum mechanical setting, we say a particle in a state ip G C£°(R n ) moves 
tangentially to F if (pV , )( ( l) G TqF for all q G F . Now, given the normal 
vectors V</? to F , we have v G T q r iff V<^(q)-v = . So a particle in a state 
ip G C^°(lR n ) moves tangentially to F when (V</? ■ pip) (x) = for all x G F , 
i.e. iff the component of the momentum of V normal to F is zero. We would 
like to generalise this notion away from C^°(lR n ) . Since F is of /U-measure 
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zero, we cannot in general specify a property for an L 2 -function ip G L 2 (M, n ) 
on r . Recall that k is a limiting procedure on the shells </? _1 [— t, t] , in which 
eventually t — ► . So it may be enough to define tangentiality on the level sets 
of ip in such a shell, and depend on k to restrict the property to Y . Fix 
a shell St '■= l P~ 1 [— to, to] around Y , then we say a ip G Dom (p\L 2 (S to )) 

has momentum tangential to the level sets of ip iff Vip ■ pip = . Since in 

dip 
terms of the local coordinate system (y, <p) of Sect. 3 this means -7— = on 

St , we see that ip \St must be constant in the normal direction (i.e. along the 

trajectories of the vector field V</? ) . So ip \St is uniquely determined by its 

restriction to Y . This defines the notion of "states with momentum tangential 

to the level sets of <p in S to ." 

Conversely, given any ip £ L 2 (Y) , we can make out of it a state i/i e L 2 (l n ) 

by extending it constantly along the normals in S to and set it equal to zero 

outside St , i.e. 

V>(x) = ip(a~^ Jx)) for x G S to and zero for x ^ S to 

(cf. proof of lemma 3.2 for a ). Denote the space of these by Tif , and the 
projection onto TiJ by Pj . Then TiJ is thought of as a "thickening" of 
L 2 (r) in L 2 (R n ) , and note that TiJ^C^ is a transverse space. However, when 
we deal with differential operators, the discontinuity at the boundary of the shell 
St i s a problem, so we prefer the following smooth version. Let Ct G C°°(R) 
be a bump function which is one on [—to, to] and zero outside [—to — s, to + e] 
for a given e . Then define for a, ip & L 2 (Y) the new thickening 

V'(x) = c i0 (^(x))-v(«; ( 1 x) (x)) . 

Denote the space spanned by these by H^ C L 2 (R n ) , then clearly Pt 7~H = 

Hi- 

We now want to constrain the Hamiltonian in such a way that it can be 

thought of as "projecting the force down to its tangential component," where a 
Hamiltonian tangential to Y will be one which keeps tangential motion to Y 
tangential, i.e. preserves the tangential states on some shell St ■ Recall that 
in the GNH-algorithm [GNH] one restricted the Hamiltonian to all states with 



-20- 

motion tangential to T . To do the same here, we can now use either of the two 
proposed methods (i) or (ii) above with the choice of transverse space as 

%o :=H C to r)C^r)DomH . 

As long as C£°(R n ) C DomH , we have that ^(% ) is dense, in which case 

this is indeed a transverse space. So assuming the latter, and that H% C C v , 
we define 

H (i) ■ K {i>) := K{Hi>) ViieT t0 . (4.1) 

In fact, when ip is smooth and H is a differential operator, this becomes 
Hu\ (ip \F) = (Hip) \F , and then the method is nothing but the one used by 
S. Helgason pp 251-252 [He] for the restriction of a differential operator to a 
submanifold. This is the first reasonable method to consider. 

To motivate the use of method (ii) , consider how a physicist might object to 
(4.1). Since H need not preserve T to , it is possible that Hip has momentum 
which is not tangential to F , in which case it seems there is a force acting on 
the particle, forcing it off F . To project this force out of the total force acting 
on the particle the restriction procedure above may not be appropriate. Instead, 
we intend to use the Hilbert space projections to first project out of H the 
parts which affect the momentum component normal to the level surfaces of ip 
in the shell </? _1 [— to, to] , then in the limit of k when the thickness of the shell 
goes to zero we will be left with the appropriate Hamiltonian on L 2 (F) (having 
maintained tangentiality of momentum to F during the limiting process). In the 
light of this we define H^ by: 

(H (ii) • k(^), «(lfe)) La(r) := Jim Y t (PfHih, hP t ih) (4-2) 

for all ipi G % for which this exists for all ip2 € % , (so that (4.2) defines a 
vector H^K(ipi) G L 2 (F) ). Since PtTt c %! ? (4-2) coincides with method 
(ii) . This is the second reasonable method. 

One can easily conceive of other methods apart from (i) and (ii) for con- 
straining dynamics, for instance forcing the particle onto F by infinite potential 
walls: 
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(iii) Given a system L 2 (R n ), <p, H as before, restrict the particle to a box 
around V , say the shell S t , but for it to be well-defined we need to as- 
sume HC^°(St) C C%°(St) and that we have a selfadjoint extension H t on 
L 2 (S t ) of the symmetric operator H \ { ip e C?(S t ) | ^(^ _1 (±t)) = } . 
In the case when H is a differential operator, this means one needs to de- 
cide how it behaves at the walls. Henceforth we assume H is a differential 
operator. If we choose ordinary reflection at the walls (quantum billiards), 
then a ip G Domi^ must satisfy 

dip 

— — = when cp = ±t , 

dip 

but without more detailed knowledge of H we do not know what additional 
boundary conditions ip must satisfy to make H t selfadjoint. So Domi^ 
is a subspace of 

V t :={il>e PtC°°(m n ) | V e L 2 (R n ), Ig = when p(x) = ±t } C <%,. 

Then in the limit when t goes to zero, the boxes force the particle (hence 
H t ) onto r . The appropriate way to take this limit, is through k , i.e. 

(H {iii) ■ k(^i), ^2)) l2{t) := }im^ (# t PWi, hP^ 2 ) (4.3) 

where ipi\St G Domi^ for all t G (0, to) and some fixed to . Note that 
the possible Vi' s are restricted by the behaviour at the boundary. So here 
for reflection, we get that (dipi/d(p)(y, t) = for all t G [— 1 , 0) U (0, t ] . 
We recognise that Pftp is in our transverse space of earlier on. When the 
space of «(V0 for such if; is dense in L 2 (T) , Hu^ is well-defined. 
Other versions of this is possible if we change the behaviour at the walls, e.g. 
introduce a phase with the reflection. 

We will not discuss any other methods, but now the surprise is that we have the 

following equivalences: 

Theorem 4.4. (1) With the assumptions above, H^ = H^ =: Hy , 

(2) moreover if H is a differential operator, then H^ = H^^ 
on DoraH( iU) . 
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Proof: (1) We first show H (i) = H {u) . Recall that C c °°(K n ) C Dom# , and 
that K(T to nC c °°(M n )) = C c °°(r) hence it suffices to show that 

(( H (i) -■ff(«))«(V'i). K (' i /'2)) L 2 (r) = 
for all tf)i E % n C c °°(M n ) . By (4.1) and (4.2), the left hand side is: 
Jim 1((I - P f T )^!, W^ 2 ) = lim 1 /" (I-P?)H1;i1>2 |Vp| dfi 

Partition T into patches U on which local coordinates exist as in 
Lemma 3.2, then the right hand side becomes a finite sum of terms: 

lim i y dp J dy(I-P?)HM2 |V^| • (*) 

(alternatively, take -i/>2 with support in U and let U vary) Now 
(I - Pj)Hi>\ =: p t J- Hj , so for all <j> e Hf D C™(S t ) we have 

= (p t , (/>)= dy<p(y) / dpp 4 

where we used the fact that on St , is independent of cp , so since 
this holds for all (which will span a dense subspace of L 2 (U) ), we 
have that J. pt(y, <p) dp = a.e. in y . Now using the continuity of 
V^IV^I , (*) becomes: 

l ™Y t j d y{^\Vp\)(y^)j dpp t = 0. 

Since this holds for all U , we find that H^ = H^ . 

(2) Next we prove that H^ = Hu^ on DomH^n) . Clearly 

P| V t \S to = Ti.f o , so it suffices to show that 
te[o,t ] 

((H {i) - HwMiln), K(ip 2 )) L2{r) =0 V^6C c °°(^)n f) DomH t 

te(o,t ] 

i.e. that lim — ((H - H t P t )tpi, hPtfa) = . 

Now ((H - H t P t )ip u hPtfa) = {(PtH - H t P t )ip u hP t ifa) , 

and let Q t =:{(/> G C?(S t ) \ <j>{p-\±t)) = } 
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which is dense in L 2 (St) , and H = Ht on Qt ■ Now given an 
open precompact set X C S t , we can write any 4> G C^°(IR n ) as 
4> — 4>x + 4>x where 4>x G Qt and fX = 0x \X . Let ^ fS^ € 
Dom H t n C c °°(S t ) for all t G (0, t ] , let G Q t and X be the 
interior of supp . Then 

{(P t H - H t P t )il;, </>) = ((P t H - H t PMx + $x), <P) 

= ((H - Httyx + (PtH - H t P t )^ x , </>) 

= ((P t H-H t P t )^ x ,4>) (t) 

Since H and H t are differential operators, they preserve supports, 
hence supp ((Pt H — H t Pt)ipx) *Z suppV'x which is disjoint from the 
support of (j) , hence since the inner product in (f) is an integral, 
we conclude that (j) is zero, i.e. (PtH — H t P t )ip _L Qt hence 
(P t H - H t P t )il; = . | 

We remark that any reasonable method of dynamics reduction should produce 

a constrained Hamiltonian H^ on L 2 (F) which is essentially selfadjoint and the 

time evolution it generates, expitH-p , must preserve the algebra of observables 

IZs when we have obtained this (see next section). We now test the above 

methods on two examples. 

Example 4.5. Recall the previous example of the sphere of radius a , F = S 2 in I 3 ; 

H = L 2 (R 3 ) = L 2 (R + , r 2 dr)®L 2 (Q) with Hamiltonian H = p 2 /2m = 

— 2~A , having domain of essential selfadjointness C£°(1R 3 ) and which 

acts on the decomposable functions V'W = f( r ) "£(#? 4>) £ C^°(1R 3 ) by 

(~^)(x) = ("^ - ? ~) f(r) ■ £(0, 0) + ^/(r) • Bm <P) 

where B = |L| 2 is the Laplace-Beltrami operator on L 2 (Q) ( O 
denotes the unit sphere). Now for t G (0, a) we have 

Hj:={iPeH | ^(x)= X[v ^ vWI] (r) •£(#,</>), £GL 2 (r)}. 

To obtain the smoothed space, choose a (^ G C£°(R) which is one on 
[—t, t] and zero on [—t — £, t + e] and set 

^:={^GH | V(x) = Ct(r 2 - a 2 ) • £(0, 0) , £ G L 2 (r) } 
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and note that if ip\T = £ G C°°(r) , then ip G C°°(R n ) C DomH . 
So we choose the transverse space % — < ip E H^ | ip\T E C°°(T) > . 
Now for method (z) for constraining H , let ip = (C,t o (p) ■ £ E % . 
Then from the explicit formula for H , we see 

where p is a function with support disjoint from </? _1 (— t, i) . So, 
since Hip is continuous near T , we have n(Hip) = (Hip) |T , i.e. 
k(H%I))(0, c/>) = (B£)(e, <p)/2ma 2 . So 

ff (i) K(V0 = «(^) = ^2 BK W 

for all ip E T t , so, since «(^t) = C°°(T) is dense in L 2 (T) , 
H({)£ = 2 ^ a2 B£, for all £ G C°°(T) is densely defined, agrees with 
the classical Hamiltonian |L| 2 /2ma 2 obtained in Sect. 2 (eq. (2.4)), 
and is essentially selfadjoint. A good result. 
Next, as an exercise in method (ii) , we calculate 

where ipi = ((t ° <p) • & as above. Now 

( ^ l)W= ( Pir 2^) f(M) 

if r G [Va 2 — t, Va 2 + 1] =: I t and zero otherwise. To work out the 
projection, we only need to consider the radial coordinate. The one 
dimensional space N generated by x 7 G L 2 (M + , r 2 dr) corresponds 
to Tif . So we need to decompose the function Xj (r)/2mr 2 according 
to N © N 1 - , i.e. i£p:Xi t ( r ) = *t • Xi t ( r ) + K r ) where A* is a 
constant and h G N 1 - , i.e. f t r 2 h(r) dr = . Now (x T , 2 mr 2 ^i ) = 
J— f T dr — At f T r 2 dr . i.e. 

2m Jit Jit 



3(v / o 2_ ^ - Va^t) 
* " 2m[(a 2 + i) 3 / 2 -(a 2 -t) 3 / 2 ] an S ° : 
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Thus {P? H^hPy-t^) 

X t • X/t (r)(S^)(^, 0) • 2n/> 2 (x)r 2 drdQ 



*Va 2 +t 

2X t I r 3 dr I (BSJ-bdn 



= 2X t ^ ((a 2 + t) 2 - (a 2 - t) 2 ) ■ (5fr, £ 2 ) L2(r) 
= 2tX t ■ (B£x, C2) L 2 (r) 

where we used the fact that the measure on T is a 2 dVt with <ifi the 
usual measure on the unit sphere. Now 

__3_r y^TT-v^ 2 "^ 1 

So * ~ 2m S (a 2 + t) 3 / 2 - (a 2 - i) 3 / 2 " 2ma 2 ' G 

Sa^( p ^^' ^ p [-t,t]^) = j™A t (Bei, 6) L2(r) 

= (^6, £2) = (#(«)&, 6) for all & G C°°(r). 

Thus -ff(ii) = B/2ma 2 = H^ = Hy on C°°(T) , in agreement with 

theorem 4.4. 
Without much extra effort, we can also show that Schrodinger operators H = 
2^-p 2 + V(q) for V smooth near T produce for £ G C°°(T) : 

H {i) Z(6, 0) = Hwffl, </>) = (w^pt) (0, (/>) + V(a, 6, 0) ■ £(0, 0). 

Example 4.6. We would like to restrict a free quantum particle in M 3 to a cylinder 
r of radius a around the z-axis . We do this in exact analogy 
with the sphere. Let the constraint be y(x) = x 2 + y 2 — a 2 , then 
V</?(x) = (2x, 2y, 0) hence the critical points are well away from F . 
So 

ip- 1 [_ Sj t] = { x G M 3 J r 2 e[a 2 - s, a 2 + t]} 

where we henceforth use cylindrical coordinates (z, r, 6) . Then 
(V>i, V^r = Jim — (-01, hP t ip 2 ) 



1 / 1^.2 



2 dt 



/ 2r 2 dr / dz / d0(^ 2 )(x) 

JWa 2 -t J-00 JO 



t=0 
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Let ipi be decomposable ipi(x) = fi(r)£i(z, 9) where & G L 2 (Y) 
and fi is continuous at r = a . Then 

(V>i,V>2)r= ft&dzM-l f ° ' r 2 Jj 2 (r)dr 

J at J^T-i t=o 

= af i (a)f2(a) / C1C2 d9dz. 

The free Hamiltonian if = —A/2m in cylindrical coordinates acts by: 
A , 10/ 0t/A 1 a 2 -!/' d 2 ip 



r dr \ dr J r 2 d6 2 dz 2 
For £ G (0, a) we have: 

Hf={^^ I iK*0 = x Jt (r) • £(*, 0) , £ e L 2 (r) } 

where It := [Va 2 — t, \J a 2 + 1] . Then with a bump function Q G 
C°°(1R) which is one on [— t, t] and zero outside [—t — s, t + s] , we 
have 

H C t :={^en J V(x) = (t(r 2 - a 2 ) • £(z, 6>) , £ G L 2 (r) } . 

Since the smooth functions of compact support are in the domain of H , 
we choose 

T t = {^eni I #rGC c °°(r)}. 

Now to obtain fl"^) , let V — (Ct ° f) ' £ £ ^~ > then 

(^)(x) = (^) Xjt (r) (L^ + ^j *(*, (?) + p(x) 



where Jt := (yaP — t, \J a 2 + t) and p has support disjoint from J t 
Now since k(Hi/j) = (Hip) |T , we get 

*««*.») = (£)(;?£ + £)«*.«) 

for all £ G C^°(r) . Thus via theorem 4.4 



2m / V a 2 09 2 dz 2 
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on C^°(r) . This is precisely the quantization one would expect of the 

constrained classical Hamiltonian (2.6). 
Remarks. (1) We regard the choice of a transverse space T as a decision on the 
direction from which H should be reduced to F . Above we chose the 
normal direction. In quantum systems with less geometry (e.g. no metric 
on an underlying manifold), it may be difficult to decide on an appropriate 
T . This choice has physical content, and seems analogous to the choice of a 
selfadjoint extension for a Hamiltonian in the Dirac approach sketched early 
this section. An alternative way of expressing the choice of the transverse 
space is to observe that on L 2 (St ) for to small enough, the pair of op- 
erators ip and P^ := id/dip form a canonical pair, i.e. [P^, 0\i/j = iip 
for all smooth ip with compact support in the interior of St ■ Then 
T\St = Ker P v fl Dom (H\L 2 (St )) ■ So in a general quantum system one 
can look for a given "local canonical conjugate" to the constraint to obtain 
a transverse space as its kernel. In fact, this observation also gives a clue on 
how to enforce second-class constraints, in the sense that if we are given a 
canonical pair cp, P v to enforce, then we do it as above, by using P v to 
locally select a transverse space to </? _1 (0) . 
(2) Recall that in Sect. 2, the classical secondary constraint we obtained for ip 
was (Vip) • p |T . The selection of the normal transverse space above, can be 
thought of as the enforcement of the constraint Wip-pip = near F , which 
looks very much like what one would expect a "quantization" of the classical 
secondary constraint to be. This leads one to ask whether we can quantize the 
infinitesimal Dirac procedure which produced this secondary constraint as a 
whole. Unfortunately this does not seem to work;- we sketch how it goes awry. 
Given the primary constraint <p and the Hamiltonian H—^p 2 + V(q) , 
the secondary constraint should be [H, <p\ "restricted" to F . Now 

[H, $[$ = -^(Av?)V - ^(Vv?) • (WO 

for all if) G C^°(lR n ) . Restriction to F is done by k . However, whilst 
— ^(V</?)•(V'^/ , ) is proportional to the desired secondary constraint Vip-ptp , 
the term —^(Aip)i/j need not vanish on F . Moreover if we use [H, (p\ 
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as a new secondary constraint to select the transverse space instead of Vy-p , 
we obtain the wrong result on the sphere. Nevertheless, given the close parallel 
which the selection of the transverse space has to the classical secondary 
constraint, we regard the use of a transverse space here as the imposition of 
a secondary quantum constraint. 

(3) We remark that one needs to reduce the dynamics infinitesimally, i.e. through 
the Hamiltonian, not directly on the time evolution unitaries exp(itH) . 

(4) In general, there is no guarantee that the constraining of a selfadjoint Hamil- 
tonian H to r produces an essentially selfadjoint operator on L 2 (F) . 
This is a difficult question that needs further investigation, and its classical 
equivalent in the GNH-algorithm also is unsolved, that is, we do not know 
whether reducing a given complete Hamiltonian vector field to a submanifold 
T produces a complete vector field on F . When F has edges, it is very 
easy to get examples of complete vector fields which will not be complete 
when constrained to F . 

(5) The confining potential approach, developed in [Ma, Fu], also treats the dy- 
namics of quantum particles on (or near) surfaces. In this, one assumes there 
is a confining potential depending only on the normal coordinates, and with a 
deep minimum on F . The system is then taken to be in a fixed eigenstate in 
the normal direction ( "normal degrees of freedom are frozen" ) , and this was 
applied to some physical examples where such confining potentials actually 
occur [Ma]. By contrast, our object is to impose the constraints exactly in 
order to understand the nature of the mathematical problems involved with 
the construction of a quantum version of Dirac's theory of constraints. As can 
be seen from this paper, such a program encounters deep pathologies which do 
not arise when the potentials only provide approximate confinement. Further- 
more if one has a physical principle which requires the particle to be exactly 
on r , or if one is defining quantum mechanics on a curved manifold, it also 
seems more appropriate to solve the constraint exactly. In addition to the 
dynamics, we also do the kinematics. 
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5. Constrained Observables. 

In this section we continue the analysis of the problem of the previous sections 
specifically in regard to the observables. That is, given the unbounded map 
K : Ctp — »■ L 2 (T) above, we wish to examine how selfadjoint operators and uni- 
taries on L 2 (R n ) lift through k to produce operators on L 2 (T) . Due to the 
unboundedness and nonclosability of k there will be some pathology, even for 
bounded operators on L 2 (IR n ) . 

The choice of field algebra JF C B(7i) will turn out to be important for 
obtaining a nontrivial constrained field algebra on L 2 (T) . In fact, the CCR 



algebra A(M 2n ) = C* { exp(iq • a) , exp(zp -a) | a G W 1 } will be too small if 
r is a curved manifold (see later in this section). A more suitable field algebra 
for curved F is Cb(IR n ) x Diff M 71 , which is here concretely the C*-algebra 
generated in B(L 2 (R n )) by the multiplication operators { T/ | / G Cb(IR n ) } , 

(7>V)(x) := /(x)^(x) , V e L 2 (R n ) = H 

and the set of unitaries { Vp I (3 G Diff lR n } , where 

(V>V)(x):=JMx) 1/2 Wx), i>eL 2 (W) 

with Jp the Jacobian of (3 G Diff W 1 . We will concentrate on these two classes 
of operators below. 

Now an operator A G B(7i) will lift through k to a densely defined 
operator A-(A) on L 2 (T) if: 

• there is a space S C C v such that AS C C^ and k(S) is dense in 

£ 2 (r) , 

• Ai\) G Ker (•, -)r = Ker k, for all tp G 5" fl Ker k . 
In this case 

A(A)K(tfj) := k(A^) VV> G 5, (5.1) 



or equivalent ly 



(A(i4)«(0), «(V0) r = Jim 1(^0, fcP t V>) VV, </> e 5. 
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Since we are interested in *-algebras of operators, we will concentrate on situations 
where the dense subspace S is invariant under the class of operators under 
consideration. Define the *-algebra: 

jr s .= | A e T | jy, e S 3 A*^ V?/> 6 S } . (5.2) 

There are three obvious dense subspaces S we can ask to be preserved, C^ , 
C c (M n ) and the transverse space T t = HJ fl C v fl domH of the last section (or 
better still, the space T^"' spanned by all % for t G (0, to] and all possible 
smoothings £ ), and we will consider all these in due course. A useful way of 
characterising T^ to ' is as 

T^ = jv>eC c (K n ) | |^(x)=o Vxe^, o<^<io}n Domif - 

The choices T t or T( to > can be thought of as enforcing secondary quantum 
constraints on the observables, cf. remark 2 of Sect. 4. However, we will see that 
in the quantum picture there is no compelling reason to do this. 

Unless we specify what S is below, we will assume some choice has been 
made. The elements of J-'s which will lift via (5.1) are those which preserve 
S n Ker/t , and for these we have: A(A)A(B)k(i(j) = K(ABip) = A{AB)k{^) 
for all ift G S , i.e. A is an algebra homomorphism. However, there are two 
pathologies associated with A ;- 
(i) given an A G B{7i) for which A(A) exists, then A(A) need not be 

bounded, 
(ii) given an A G B(J~i) for which both A(A) and A(A*) exist, we need not 

have that A(A*) C A(A)* . 
5.3. Example of (i) : 

Let 7i = L 2 (R 2 ) , y?(x) = X2 so V = </? _1 (0) is the xi-axis. Choose 
S = C c (R n ) . Clearly all V p preserve S = C C (R 2 ) for all (3 G Diff R 2 , so 
since (V^V)(x) = J' 1 ' 2 ^' 1 ^) ■ V>0# _1 x) we have V p G T s • If £T C T , we 
get that Vg preserves C c (M 2 \r) = Ker (•, -) r n C C (R 2 ) , and so A(Vp) exists. 
Consider now the (3 G Diff R 2 given by fl(xi, x^) '■= (axi, X2e Xl ) , a > , so 
clearly /3r = T and Jg(x) = ae Xl 7^ . (That (3 is a diffeomorphism is clear 
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since its inverse is f3~ 1 (x\, X2) = (x\/a, X2e~ Xl ' a ) which is also differentiable) . 

Now for f G C C (I 2 ) we have 

(A(V^)k(VO, A(V»^)) L2(r) = (V^, V^) T = (V^\T, V^\Y) L2{T) 

|(Wrr)| 2 dxi 

Jp(xi, Q)\ip(axi, 0)| rfxi 
- / e^^x, 0)| 2 dx=: .fy 



and now we can choose a ?/> G C C (IR 2 ) with ||«(V , )|| = 1 1 V 7 1^11^2^ = 1 which 
can make 1^ arbitrary large, e.g. V — X\ n n +n x [o 11 anc ^ cons ider 1^ as 
n — ► 00 . Thus A(Vg) is unbounded. 

We remark that if we chose a more restrictive space £ , e.g. T^ to ' , then 
the diffeomorphism in the last example will preserve T^* -* , so this pathology 
cannot be removed by enforcing "secondary quantum constraints." 
Example of (u) : 

Continue the previous example, noting that A(V^) exists because (3~ 1 T = F . 
However, if A(V£) C A(V£)* we have for all ip G C C (K 2 ) that 

(A(^)k(VO, A(V p )k^)) = (k(V0, A(^)*k(W)) 

l (ij \ 'L 2 (r) 

= (k(V), A(y;)«(^)) La(r) 
= («(^), «(y;^v)) La(r) 

= («(^), «(V0) L 2 (r) , 

which makes A(V^) unitary, hence bounded, in contradiction with example 
{%) , so it is false that A(V£) C A(V r /3 )* . 

Now, it appears reasonable to the authors that in a constraining method, 
boundedness and the adjoint operation should be preserved, at least on the physical 
variables. So we want to restrict the set of operators under consideration to those 
satisfying these two requirements. Observe that an A G J-'s will lift to a bounded 
operator A (A) on L 2 (T) iff 

(Aip, Aip) r < M ■ (V>, V)r Vi^g5 and a fixed M < 00 . 
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(Clearly if ip G S fl Ker k then the last inequality implies (Aip, Aip)r = , i.e. 

Aip G Ker k , so A (A) exists). 

Definition: 

O s := {Aef s \ (AtP, AiP) r < M A • (</>, V)r > (^*V>, ^*V0r 

V^6 5 and some M A < oo , and A(A*) C A(A)* } 
D s :={ie0 s | AS CKer kDA*S} . 

Note that if A G O5 , then A (A) is bounded, so we have in fact that A (A) = 
A (A*) , so there are no problems with selfadjointness if A is selfadjoint. These 
will occur however if we consider A on unbounded operators. 

Lemma 5.4. Os and T>s are *-algebras. 

Proof: Both sets are linear spaces. To see that Os is closed under taking of 
adjoints, take the adjoint of A(A*) C A(A)* to get A(A**) = A(A) C 



A(A)=A(A) CA(A*) , i.e. A* is in C s if A is. Then clearly 
T>s is also closed under taking of adjoints. To see that Os is an 
algebra, let A, B G Os , then 

(ABip, ABip) T < M A ■ (Bip, Bip) v < M A M B (ip, ip)v 

for all ip G S , and similarly for B* A* . Moreover 

A(AB)* = [A(A)A(B)}* D A(B)* A(A)* D A(B*)A(A*) = A((AB)*) 

and thus AB G Os ■ That T>s is an algebra is obvious. | 

Now V s = Ker (A\O s ) is a *-ideal of £> s , and TZ S := 0s/£>s = A(C S ) C 
i3(L 2 (r)) . We think of 7£s as the "physical observables" obtained from enforc- 
ing the constraint on Os , in analogy to the algebra TZ of the T-procedure for 
Dirac constraining. IZs cannot be zero because the identity operator I G Os 



and A(l) = I . We do not expect IZs to be a C*-algebra, but it can easily 
generate a C*-algebra since 1Z S = A(O s ) C B(L 2 (T)) . 

Given that the commutant ip' is the traditional observables, we find here: 

Theorem 5.5. Fs n 0' Dti C O s ■ 
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Proof: Let A G J-'s H (p' C\ h! , which is a *-algebra, so it also contains A* 
Then 

(Aip, hP t Aif>) = Wh^PtA^f = \\AwT 2 P t ipf 

< \\Af ■ \\hV 2 P t ^\\ 2 = \\A\\ 2 ■ (ip, hP t 1>) 

for all ipeS . Thus 

(Aip, Aip) T = lim Ji(Ai(), hP t Aip) 

< || A|| 2 lim ^(^, hP t *P) = || A|| 2 • (</>, V)r • 

The same is true for A* . Now for all ip, <p E S : 



(A(A)k(VO, «(0) La(r) = (^W)' "(M^r) 
= (A^, 0) r = lim^W, ^*0) 

= lim ±M), A*hP t (/>) = lim ±-M, hP t A*(j>) 
= (V, A» r = (k(V0, A(A*)«(0)) L2(r) 
so A(A*) C A (A)* . Thus A e O s . | 

Remarks. (1) Note that due to the limit in (•, -)p , it is only the behaviour near Y 
which contributes in the constraining of an operator. In fact the proof actually 
shows that an A e J-'s will be in Os if 

]hn±\\h*? 2 P t Aib\\ 2 = lim UAhV~ 2 P t ib\\ 2 



and lim±(V>, [A*, kP t ]0) =0 

for all ip, (p E S and the same for A* . So if A commutes with (p and 
h on a small neighbourhood of Y , it will also be in Os ■ Behaviour of A 
away from Y is irrelevant. 
(2) We note by 5.5 that all multiplication operators by bounded Borel functions 
which preserve S will be in Os ■ In particular, for a reasonable physical 
system we expect (p to preserve S , so the constraint (p is in Os ■ In 
this case (p is in T>s by the following argument: 

lim||/^ 2 P t ^|| 2 = limi||<^V 2 PW'l| 2 

<limisup{|v?(x)| 2 | xev^M-MH^iW 

= (limsup{|^(x)| 2 | xe^K t]}).limi||/^ P ^||2 

=0 V^eS 
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using the facts that the limit lim j||/i 1 / 2 P i '0|| 2 = lim.j(ip, hPfip) exists 
due to ip G S C C v , and that the limit of the supremum is zero. Hence 
((pip, (pip)r = , i.e. (pip G Ker k for all ip G S , so A(<£) = as 
expected. 

Theorem 5.6. An A G Ts preserves Ker k fl S iff 

limJ\\P t A(I-P t )iP\\ 2 = 

for all ip G Ker kC\ S . 

Proof: We need the following lemma: 

Lemma: With assumptions and notation as above, we have 

ipeKer(-,-) r iff lim j\\P t ip\\ 2 = . 

Proof: ip G Ker (•, -) r iff = (ip, ip) T = lim ±(ip, hP t ip) . Now 

(ip,hP t ip)= f |Vv?(x)HV(x)| 2 dMx) 

>inf{|V^(x)| | xe^I-M]} / |V(x)| 2 ^(x). 

Jvp-i[-t,t] 

Since |Vy?(x)| is assumed to be bounded, continuous and nonzero on a 
neighbourhood of F we have 

liminf { |V^(x)| | x G ^ _1 [-t, t] } =: M > . 

1 If 

Thus lim -(?/>, /iP t V) > M lim - / |^(x)| 2 rf^(x) 



^°t ^°t j<p-i[- t ,t] 

= M lim -\\P t ip\\ 2 

so ip E K.ei (■ , -)r implies = lim —(ip, hPtip) > M lim — \\Ptip\\ 

L r \J (y L '\J ly 

and so: lim — ||P+'0|| = 

On the other hand, given a ip G 7i satisfying the last equation, we have 

(ip,hp t ip)= [ |v^x)HiKx)| 2 d/i(x) 

Jip-^[-t,t\ 

< sup{ |Vv?(x)| | x G (p-'l-t, t] } [ |^(x)| 2 ^(x) 
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and the the limit of the supremum (denoted N ) exists by the bound- 
edness assumption, and is strictly positive. So 

1 1 

lim -{ib, hP t ib) < N ■ lim -\\P t ip\\ 2 = 

i.e. lim j (ip, hPtip) = , i.e. ^eKer(-,-)r . ▼ 



So by this lemma, A G jFg preserves Ker kDS iff lim t ||PtA-0|| 2 = 
for all tb G Ker kC\ S . Now 



lim UP t AP t i/>\\ 2 < lim ^PtAPtfWP 



t"\\ 



< ||A|| 2 -lim±||P^|| 2 = 



t^o 



whenever ib G Ker k . So by the triangle inequality: 

\\P t AP t iJj\\ + \\P t A(I - P t )ip\\ > \\P t AiP\\ 

> \\\P t AP t iP\\-\\P t A(I-P t) 
we obtain: lim 4-||P t ^|| = lim 4-||P t A(I - P t )ib\\ 

for all ip G Ker k(~] S . | 

An immediate consequence is that A(A) exists for all A G (p 1 fl J-'s (though 
A(A) need not be bounded in general), and as before, we only need commu- 
tativity close to T to get this. Also note the similarity with the condition 
Pph ys ^4(I — -Pphys) — for an observable in Dirac constraining. Moreover the 
proof of 5.6 used only the existence of the limit (3.6) for ib , so it will work for 
larger domains than C v . 

In the particular cases where £ is either C v or C c (lR n ) , we shorten the 
notation to A<p := A c and A c := A c (Rn ^ where A can be J 7 , O, V or 
TZ . Note that 

C c (R n ) n Ker k = { $ G C c (R n ) \ $ \Y = } = C c (M n \r) 

V c = { A G O c | AC c (R n ) C C c (M n \r) D A*C c (R n ) } 

and for an A G O c we have A(A)k(iJ)) = K,(Aip) = (Aip)\T for all ip G 
C c (R n ) . Analogous statements hold if we replace C c (lR n ) by C v . 
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On comparing the method above to that of Landsman [La], we note that the 
algebra which Landsman selects to impose the constraint on is the subalgebra of 
the commutant ip' which preserves C c (W n ) , i.e. ip' fl T c . In our case the 
algebra which we constrain, O c can be considerably larger than that, with con- 
sequently larger algebra of observables TZ C , but on the other hand there are also 
nonzero A G <p' fl T<\O c (see below). 

Next we wish to examine whether particular classes of operators are in Os ■ 

• Consider the multiplication operator (Tfip) (x) := /(x) ■ V>(x) \/ip E Ti 
where / is bounded and Borel. It is not automatic that Tf G J-'s • In 
fact if S = C c (R n ) , then T f C c (R n ) C C c (R n ) iff / is continuous (in 
which case Tf G O c by 5.5). For the choice S = C v we see that if a 
Borel function / is discontinuous on F , then restriction to F may not 
be defined, i.e. we may have Tf G" T^ . (This also shows 'p'\J r ip 7^ hence 
that Lp' (jL 0^ U O c ). On the other hand, if / is continuous and bounded 
on some shell St , we have Tf G T^ . So T v contains a larger class of 
multiplication operators than JF C , and these are all in O^ . In particular, 
if exp iq ■ a G T then exp iq • a G O c fl O^ . However a Tf is only in 
T r w if / e T« . 

• Now consider the unitaries Vp , f3 G Diff IR n . For the choice S = C c (M. n ) , 
we have that both Vp and VS preserve C c (lR n ) , hence Vp G T c for all 
[3 G Diff W 1 . This is not however true for the choice S = C^ . In fact, let 
[3 be a fixed translation /3x = x + a not preserving F , and let ip G C^ 
be continuous on a neighbourhood of F but so discontinuous on (3F that 
restriction to (3F is not defined. Then Vp-iifj G" C^ , and so Vp-i G" T^ . 
Thus T^ contains a smaller set of the unitaries Vp than T c . Below we 
will examine when these are in O v . 

• Let S = C c (R n ) 3 ip , then 

A(Vp) k(v) = «(W) = (W) tr = (J^ /2 rr) • (^ o /?) rr . 



So A(^)=I iff j8x = x Vx G T and Jp\F = 1 , i.e. Vp-IeV c iff 
/3x = x Vx G T and Jp\F = 1 . 
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Since k(VO = V'lT for -0 G C c (IR n ) , one might surmise that VpY C r iff 
Vp G C c . However, the example (5.3) shows that it is not true that j3Y CT4 
Vp E O c . For the converse, we have: 

Lemma 5.7. If VpeO c l)0 9 , then (3Y C T . 

Proof: Let V£ G C c , hence Vp preserves Ker (•, -)r n C c (R n ) = C c (M n \r) . 
This is equivalent to /3(R n \r) C R n \Y , which is equivalent to (3Y C 
T . This argument directly adapts to C^ . | 

Note that if Y is curved and nonperiodic, then the translations will not preserve 



it, so exp(zp -a) G" O c U O v hence of the generating unitaries of A(M 2n ) , only 
the commutative set { exp(iq -a) I a G M n } is in O c U O v . This is why we 
consider the CCR-algebra as too small a choice of field algebra for this type of 
constraining. 

Theorem 5.8. Given notation above, we have: 

VpeOcUO^ iff [3YCY and J^(x) = J^(x) for all xGT 
where Jp (resp. Sp ) is the Jacobian of (3 (resp. (3\Y ) 
with respect to \i (resp. 7 ). 



Proof: Let Vp G O c U O v , so by 5.7 (3Y C Y . Now since A(V$) = A(Vp)" 
we get 

A(V>)*A(V>)k(V>) = A(V>)*k(W) 

= A(Vp*)K{Vp^) = K{Vp*Vpi}) 

= K{il}) = A(Vp)A(Vp)*K{il)) V^GC C (R") or C 9 . 



Thus A(Vp) is unitary, so 

(V>, V)r = (VpiP, Vp^) T = (Vp*P\Y, Vpij\Y) L2(T) 

= J |(X 7J V)(x)| 2 rf 7 (x) = J J^x) • |V(/3x)| 2 d 7 (x) 

= 1 J^(x)-J j ; 1 (x)-|^(/3x)| 2 d7(/3x) 

= J Jpip- 1 *)-^ 1 ^- 1 *) ■ |^(x)| 2 rf 7 (x) 

= /"|^(x)| 2 rf 7 (x) VVGC C (K") or C^ 
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iff Jpifi- 1 *) ■ Jp 1 ^- 1 *.) = 1 for all x G T , i.e. J^x) = JJ^x) for 
all x G T . 

Conversely, observe that (3T C T guarantees that both A(Vg) and 
A(Vg) exist, so by the reversibility of the previous calculation when 
Jp — J/3 on T , we see that A(Vg) is unitary, hence bounded, and 
similarly the same is true for A(V^) . Since A is a homomorphism and 
inverses are unique, we get that A(V^) = A(Vp) , i.e. Vp G O c UO v . 



Remarks: (1) Note that if (3 G Diff M. n preserves all the level sets of ip we have 
ip o (3 = ip , so [V/3, <p] = . However since we need not have Jp = Jp on 
T , there are certainly such (3 for which Vp £ O c . Thus jF c n^'\(9 c ^ , 
i.e. Landsman [La] quantizes some operators which we have excluded from 
our observables. (Some of these can still be taken through A using 5.6, 
but they need not preserve the adjoint or boundedness) . Nevertheless, we 
conclude that the field algebra Cb(M. n ) x Diff M. n has ample elements in its 
physical algebras 1Z C or IZ^ . 

(2) There is no reason in general to expect that the time evolutions will preserve 
IZs , so we will need to extend it by taking the constrained field algebra as the 
C*-algebra generated in £>(L 2 (T)) by IZs an d the unitaries expitHr , 
fGl. 

(3) Recall that we consider the enforcement of secondary quantum constraints as 
the selection of the set S = T^ to ' which the observables should preserve. 
Since this is an analogy to the classical procedure, it is natural to look for 
reasons to justify such a choice, and a first attempt may be to ask whether 
the choice of transverse states for S will make the Jacobian condition in 
5.8 obsolete, i.e. whether for a Vp G T T (t) we have Vp G T (t) iff (3T C 
r . This is not true, as we can see by an easy counterexample. Continue 
the example (5.3) with the additional assumption that the Hamiltonian has 
DoiaH = C^°(1R 2 ) and that (3 G Diff M 2 is not the given one in (5.3), but 
/3(xi, X2) = (%i, 2^2) • Then both Vp and VS preserve 

T(*°> = jveCf(R n ) I J^( x )=° Vxg^, t^<to\ 
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hence Vp G J r T m ) ) , and j3Y — V , /?|T = I , Jp = 1 , Jp = 2 , hence 

so Vp £ O T (t 0) . 

Another possible reason one may want to use to justify the enforcement of 
secondary quantum constraints, is to produce a common dense invariant do- 
main k(S) for the observables IZs • However, for the choice S = T^> 
there is no reason why the Hamiltonian Hr should preserve k{T^">) , so 
it will be unreasonable to require the observables to do so. 
So at this stage we fail to see why secondary quantum constraints should be 
imposed on the observables. 
(4) It is interesting to observe that for any Hilbert space operator K : Tii — > 
7^2 with dense range, we can obtain the structure above. That is, given a 
dense space S C domif for which K(S) is dense in H2 , we can define 
jFg, Os, T>s, IZs exactly as before by replacing T with B(7ii) , and k 
by K . So every such operator K and space S defines a lifting problem, 
hence a short exact sequence of *-algebras 

-► v s -► o s -+ n s -* . 

This then produces a similar short exact sequence for the automorphism 
groups and their lifting through the factorisation. In particular, for ordinary 
Dirac constraining the operator K is the projection P p h ys : 7~t — ► 7i p hy S , 
the space S is TL and As^J~s is just A where A denotes either O 
or T> . This way of thinking nicely unifies the current construction (where 
K — k ) with that of the Dirac approach. The main difference is that -Pphys 
is bounded whilst k is unbounded and nonclosable. 
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6. Constraining by a general selfadjoint operator. 

For later reference we start by summarizing the constraining algorithm developed 
in the preceding sections. Given the following data: operators p , q and H 
on C c °°(M n ) C L 2 (R n ) , a bounded <p e C°°(R n ) such that T = </? _1 (0) is a 
smooth (n — 1) -dimensional submanifold and moreover V(p is bounded and 
nonzero on a neighbourhood of F and a unital field algebra JF C B(L 2 (M. n )) 
containing (p , exp ip • a , exp iq ■ a , 

(1) define (V>i, ^)r := Jim i (V>i, hPtfo) W>* G C^R"), 

L 2 (r) = [C ¥ ,(R n )/Ker(-,-)r] — (v)r , k : C^ ^ L 2 (T) is the factorisation 
map and L 2 (F) is the constrained space. 

(2) Let T to :=Ti.5r\C<priDomH , check that HT to cC^ for t small enough, 
and if so, define the constrained Hamiltonian Hr by H-p ■ k(V0 := k(Hi/j) 
for all if) G 7t . 

(3) Choose a space 5 in 7i for which k(S) is dense (e.g. C v ), and define 

J r s = {AeJ r | ASCSDA*S} 

O s = {Aefs | W, M) T < M A (^, *p) r > (A*V, A*^)r, A(A*) C A(A)* } 

where A.(A) denotes the closure of the operator defined by 

A(A)«(^) = K{Ai>) V^e5. 

Define V s :={AeO s | ASCKer^D A*S } , then 

A(Os) — Os/T>s =: Us 5 an d the algebra of the constrained observables on 
L 2 (r) is the C*-algebr a generated by IZs and exp(zRifr) ■ 
In this section we want to generalise the method above to impose a constraint 
" Cip — " where C is a general bounded selfadjoint operator on a Hilbert 
space 7i with zero in its continuous spectrum. (If we start with an unbounded 
selfadjoint operator C we can, without loss of generality convert the problem to 
a bounded one by replacing C with f(C) where / is a continuous bounded 
real-valued function with f(x) = x on a neighbourhood of zero). As before, 
we still assume that there is a unital field algebra T acting on Ti , containing 
C , and that there is also a (possibly unbounded) Hamiltonian H given on 
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7i , ex.p(itH) E T for all t . We will be concerned with the construction of 
three objects;- the constrained Hilbert space (and the constraining map to it from 
the original space), the constrained Hamiltonian, and the algebra of constrained 
observables. 

At the abstract C*-level, the kinematics part has already been solved [GH], 
but one obtained the set of all representations in which the constraint can be 
imposed as an eigenvalue condition. Here we want to build a particular concrete 
constrained system out of the given unconstrained one. 

Now geometry was paramount in the analysis of the previous sections, in that 
we needed a metric (on an underlying space) to define gradients, norms of vectors, 
normals to surfaces and the Lebesgue measure. This vital piece of information is 
missing in the problem under consideration, and we somehow need to augment 
the given data {JF, 7i, H, C} in order to adapt the method previously found to 
this problem. The extra information we will assume is: 
(i) a maximally commutative C*-algebra A C B(7i) containing C . 

(Call A a polarisation, and we know that it always has a cyclic and sepa- 
rating vector [BR 2.5.3]) 
(ii) A choice of a cyclic and separating vector Q for A ■ (Call O the vacuum), 
(iii) A choice of scaling operator K E A+ , Ker K = {0} . 
(iv) A selfadjoint operator Vc such that [Pc, C]ip = iip for all ip E V = a 
dense invariant subspace of P t H where P t denotes the spectral projection 
of C of the interval [—t, t] . This is thought of as a "local" canonical 
momentum for C , which will define the normal direction to the constrained 
system. 
Now define the space C := { A E A | ujo(A) exists } where 

(PA AP t il) 

Since obviously I and P t E £ , this space is not zero. Moreover, the space 
is selfadjoint, and we also have that |u;o(^4)| < ||A|| for all A E £ because 
|(Pt^, AP t Q) j || P t Vt || 2 1 < ||A|| , and so ujq extends to the closure of £ . Now 
since u$ is obviously positive on the positive elements of £ , we can apply [KR] 
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4.3.13 to conclude that ujq extends to a state on A . Henceforth we fix a choice 
of extension and still denote it by ujq . 

Denote the GNS-representation of ujq by (tto, Oo, T~Lq) ■ Noting that 
ujq{K) > , we define a state uk on A by 

u K (A) := u (KA)/u (K) , AeA. 

Then, inspired by the following lemma and subsequent example, we identify the 
constrained Hilbert space with Ti^ 

Lemma 6.1. lok is a "Dirac state" in the sense that ujk(AC) = for all 
AeA, i.e. n (C)Q = . 

UJ K V ' UJ K 

Proof: It suffices to show that uk(C 2 ) = , since C is selfadjoint. So 

u (KC 2 ) _ (^ KP t cri) 
u (K) \\P t npu (K) ' 

Now \(KP t n,p t c 2 n)\ < \\kc\\ ■ \\p t n\\ ■ \\p t cn\\ 

= \\kc\\- \\p t n\\ ■ || / \dP(\)n\\ 



-t 

< \\KC\\ • ||P t fi|| -t\\ I dP{X)Q.\\ 



= t\\KC\\ ■ WPfVt 



-t 

2 



Hence \uJ K (C 2 )\ < — k* limt||KC|| = . I 

Remark: The proof above is easily adapted to show that uq is also a Dirac state 
on A . More precisely, we have for the left kernels that iV Wo = N L0K due 
to Ker K = {0} and the fact that A is commutative. Hence we have an 
identification map 9 : A/ N UJo — ► A/ N LJK by 6>(A + A^ ) = A + iV WK and 
6* extends to a map 6* : 7i WQ — ► 7Y WK because 



1/2 

w K (i*A) 1/2 = wo(^A) 1/2 



||0(A + N W0 )\\ = u K {9(A + N U0 )* 6(A + N Uo )) 



< \\K\\ -uo(A*A) 1/2 = \\K\\ ■ \\A + N ( 



woII-Wq- 



To define a constraining map k : 7i — *■ Ti from the unconstrained to the 



UJK 



constrained space, recall that O separates A , hence the map 

k(AQ) := ir UK (A)Q UK , AeA 
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is well-defined on the dense subspace Afl with dense range. This will be our 

choice of constraining map for this context. Since ||k(AO)|| 2 = ujk(A*A) , we 

see that Ker k = N LUK O = N Uo O . 

Example. To motivate the preceding structures, and to see what is involved in the choices 

of A, fi, K, we now reconsider the constraint situation ip , L 2 (IR n ) of 

the preceding sections. 

Starting with ip = (p(q) , the natural choice of polarisation is A : = 
{/(q) | / G L°°(i") } = (exp(zq-a) | a£l n }" , and for O we can 
then choose any positive nowhere vanishing L 2 -function . Choose the Gaus- 
sian O(x) = exp(-a|x| 2 ) , a > and note that C c (R n ) = C c (R n )0 , 
hence C c (IR n ) C AVt . Now for K we must have K G A+ , so K — fc(q) 
for some k G L^(IR n ) , /c _1 (0) = . For the moment we will choose k 
to be continuous, and then below deduce the precise choice which will cor- 
respond with the previous results for this situation. In particular, what we 
want to show is that 

( K (ACl) t k(BCI)) = b((AQ\T), (Snm) (*) 

n ulK l (i ) 

for all A, B G C c (M. n ) = C for the right choice of K , where k is the 
map defined in this section and b is a normalising constant. Now 

(Aft, KP t BQ) 
(.(An), «<«!)) ^ = UK( A>B) = Km \ mpMK > 

whenever A*Be£ (which we will see below is all C c (IR n ) ) so if we write 
A = //l(q), B = /s(q) with j a-, Ib G C c (R n ) and remember that 

(PtV)(x) =Xp-i[_ M ](x)Xx) , 

we have 

(AQ, KP t BQ) ^_ 1[ _ ti] 7^(x)A;(x)/ B (x) exp(-2a|x| 2 ) dfi(x) 



\\ P M\ 2 J^-i[_ M] exp(-2a|x| 2 ) d/z(x) 

Now recall that for / G C c (IR n ) and small t we have 

/ /(x) <fc(x) = T ( f f { * ,(P \, <h v (yj) dp> 
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J |V(y, 0)| 
= 2tJ (f/\V<p\)\Yd 1 + 0(t). 



i / s , ^ (Afi, KP t BQ) 

thus: (.(AQ), «(fln))^ = Jim \ m?uJo{K l 

jJlAkfBe/lV^lTd^ 



oj (K)J r (e/\V ( p\)\Td 1 

where e denotes the function exp(— 2a|x| 2 ) . Now since by a similar argu- 
ment 

(p,n, A-p,n) J r (fa/|v y |)trd 7 
<*(*> = >™ p^p = /r(e/ |v^ r r*y 

(the existence of this limit shows C D C c (M. n ) ) we have 

f r (7 A kf B e/\V(p\)\Yd 1 



(k(ASI), k(BQ)) 



n„ K j T {ke/\V^\)\Td 1 



iA7 A fBe)\rd 1 ((Am,(Bm)_ 



L"(T) 



f r e\Td 7 \\n\r\\ 2 L2{r) 

(where we made the choice k = \S/ip\ ), which is just the normalised inner 

product of L 2 (Y) , and this will produce the same k than we had before. 

If there is also a group of symmetries which needs to be unitarily implemented on 

the constrained Hilbert space 7i LUK , we can use this unitarity to determine K . 

So now that we have the constraining map k : 7i — ► 7i as above, we 

can proceed to constrain the dynamics. This is where we will use the assumed 

operator Vc in the role of the "normal derivative i-B- near Y ." Thus start 

by selecting a space of transverse states by 



Hj:={P t ij I VeDomPc, P t V c ^ = 0} 
and % := { ip 6 DomH \ P t if; eUj) 

and note that % C T s if t > s . For a consistent constraining we need to 
assume there is a to > such that the set St := {ip £ % | Hip E .4.0} 
satisfies Ker k n St () = {0} and K>(St ) is dense in 7i LUK . Given this, we 
define the constrained Hamiltonian He on domain K(St ) to be 

H c ■ k(V) == <Hijj) V^e5 t0 . 
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Example. Continue the previous example with C — ip and choice K = \S/ip\ . 
Choose 

V^:=if{ip)^ ty>eCHR n ) 

where / is a smooth bump function which is one on [—to, to] and 
zero outside [—to — e, to + e] for some e > . Then 

H — {ip G H | suppV' C <£> _1 ([— to, to]) , ip(y,<p) independent oft/? 

on <£ -1 [-t ,t ] } 
If H is a differential operator with domain C^°(M. n ) , we have 
T to = { V> e C^°(IR n ) | ^[^ _1 ((-t , t )) is constant in the direction V<p } . 
Since U\ n C c °°(M n ) C % , we conclude that 

H^Kty) = K(Hif>) WeT to 

defines the same operator as before. 
Finally, to obtain the constrained kinematics as in Sect. 4, let S C Afl be 
a subspace such that k(S) is dense in 7i U}K . Define in complete analogy with 
Sect. 4, 

Ts := { A e T | Aip e S 3 A*^ V V e 5 } 

O s := { A g ^s | (A^, A^)k < M A • (V>, Vk > (A>, A» K W> g 5, 

M A < oo and A(A*) C A(A)* } 
V s :={AeO s | AS CKer kDA*S} 
where A(A) k(V) := /c(#) Vt/> G 5, A G C s • 

The previous proofs carry over verbatim, so T>s — Ker (A\Os) and the algebra 
of observables acting on 7i U}K is the C*-algebra generated by expiRHc and 



U s := Os/Vs = A(O s ) C B{H 



UK . 



Example. In continuation of the last example, it is obvious that if we take T = 
C b (R n ) x Diff R n and S = C c (R n ) , then we will get the same observ- 
able algebra than before. 
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Remark Another way to proceed could have been to extend the Dirac state ujk 
on A to T . However, this would have enlarged the GNS-space (pro- 
ducing a wrong result for the example of this paper), and in addition 
there is the problem that the extension need not be unique. 

7. Conclusions and Discussion. 

To summarize;- we took the problem of restricting a quantum particle in M. n to 
a lower dimensional submanifold not conserved by the dynamics. We found a 
method to do the restriction, and this came in three parts;- a map k from a 
dense subspace of L 2 (R n ) to the constrained Hilbert space L 2 (T) , a method 
for constructing the constrained Hamiltonian from the unconstrained one, and a 
method for the construction of the algebra of observables TZr on L 2 (F) . For 
each of these three parts there were choices to be made;- for the first part a domain 
for k ( C c (lR n ) and C v seem the canonical choices), for the second part a 
choice of transverse space (here tangentiality of momentum to the level sets of ip 
near F seems to be the physically justified criterion), and for the third part a 
choice of space S on which to reduce the observables (likely choices here are 
C c (M. n ) , C,p or the transverse space used) . We regard these choices to be made 
as physical ones, and in the two examples, the choices made were evident, and 
produced results in agreement with the known solutions. The necessity of these 
choices should be compared with the choice of boundary conditions required in 
the constraining of a free quantum particle to a box. 

There are mathematical pathologies associated with our proposed method, of 
which the central one is that the map k is neither bounded nor even closable. As 
a consequence, the constraining of observables need not preserve boundedness or 
involution, so we had to restrict the method to those observables for which these 
pathologies do not occur. The method is also sensitive with respect to the original 
choice of field algebra, and may not provide enough observables for some choices. 

Whilst this method was developed in analogy with the classical constraint 
method, the analogy became quite vague at certain points, for instance, one can 
interpret the neccessity for a choice of transverse space as an imposition of a 
secondary constraint, but this constraint could not be interpreted as obtained 
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from a similar method as the classical one. Moreover secondary constraints do not 
seem to be neccessary for the constraining of the observables. 

We expect from our work in [GH] that the pathologies which occurred here 
can be circumvented by a suitable generalisation to a C*-algebra framework, inso- 
far the T-procedure solves already without pathology the kinematics part of the 
question there, but this would involve losing the original representation which may 
contain some physical information. The advantage of the current method is that it 
constructs the constrained system concretely from the given representation of the 
original system. Nevertheless, we intend to further pursue this line of thought. 
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